The high-pressure superconductivity in SiH 4 : the strong-coupling approach 
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In the paper, the thermodynamic parameters of the high-pressure superconducting state in 
the SiEU compound have been determined (p = 250 GPa). By using the Eliashberg equations 
in the mixed representation, the critical temperature, the energy gap, and the electron effective 
mass have been calculated. It has been stated that the critical temperature (Tc) decreases from 
51.65 K to 20.62 K, if the Coulomb pseudopotential increases (y* 6 (0.1,0.3)). The dimensionless 
ratio 2A (0) /fcsTc decreases from 4.10 to 3.84, where the symbol A (0) denotes the value 
of the order parameter close to the zero temperature. The ratio of the electron effective mass 
to the band electron mass is high, and it reaches maximum equal to 1.95 for the critical temperature. 
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The metallic state in the hydrogen is intensely studied 
since 1935 pj. It has been shown that the hydrogen is 
metallic under the high compression (~ 400 GPa) [2[. 

At the low value of the temperature, the metallic 
state of the hydrogen transforms into the superconduct- 
ing state In the pressure (p) range from ~ 400 GPa 
to ~ 500 GPa, the value of the critical temperature (Tc) 
reaches maximum equal to 242 K (p = 450 GPa) |J, Q . 
For the pressure's values from ~ 500 GPa to ~ 800 GPa, 
the critical temperature increases: Tc € (282, 360) K 
0, Q . The extremely high value of the critical temper- 
ature has been predicted near 2000 GPa, where Tc can 
assumes even ~ 600 K 0, Q . 

Recently, it has been suggested that the hydrogen- 
rich compounds become metallic at the pressure's value 
lower than for the pour hydrogen. In particular, the fol- 
lowing compounds have been taken under consideration: 
methane (CH 4 ) fl(|, silane (SiH 4 , SiH 4 (H 2 ) 2 ) QjHll, 
disilane (Si2H 6 ) [14j . and germane (GeH 4 , GeH 4 (H 2 )2) 

mm. 

In the case of SiH 4 (H 2 )2 under pressure at 270 GPa 
and Si2H6 at p = 250 GPa, the theoretical results have 
predicted the existence of the high-temperature super- 
conducting state: Jh = 129.83 K and T c = 173.36 K, 
respectively [l7[, [3. We notice that the above values 
of the critical temperature are even higher than for the 
cuprates fioj ] . 

In the presented paper, we have studied the thermo- 
dynamic properties of the superconducting state in the 
SiH 4 compound under the pressure at 250 GPa. 

We notice that the experimental data have confirmed 
the metallization in SiH 4 (p ~ 50 GPa) [HEl. Further- 
more, the critical temperature increases with the pressure 
and assumes the maximum equal to 17.5 K at 96 GPa and 
at 120 GPa. 

In the SiH 4 compound (p = 250 GPa), the strong 
electron-phonon interaction has been predicted [2(|. In 
particular, the electron-phonon coupling constant takes 



the value: A = 0.91. Taking into consideration, the above 
result one should expect the high value of the critical 
temperature. In the considered case, the remaining ther- 
modynamic parameters are probably beyond the BCS 
predictions [Hi, HI- For this reason, the numerical cal- 
culations have been made in the framework of the Eliash- 
berg approach (23|. 

The Eliashberg equations defined both on the real and 
imaginary axis (the mixed representation) can be written 
in the following form [24| , (2|| : 
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The symbols <p (w) and Z (uj) denote the order parame- 
ter function and the wave function renormalization factor 
on the real axis (uj); <p, n = <fi(iuj m ) and Z m = Z (iuj rn ) 
represent the values of these functions on the imaginary 
axis (i = v— I). The Matsubara frequency is given by: 
ui m = (n/f3) (2m — 1), where /3 = (ksT)~ (ks is the 
Boltzmann constant). Both on the real and imaginary 
axis, the order parameter is defined as: A = <j>/Z. 

The electron-phonon pairing kernel has the form: 
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where the Eliashberg function (a 2 F (f2)) has been calcu- 
lated in the paper [20(; the maximum phonon frequency 
(!l max ) is equal to 426.1 meV. 

The function fj,* (u> m ) = /i*9 (ui c — \uj m \) describes the 
electron depairing interaction; /Lt* is the Coulomb pseu- 
dopotcntial. Due to the absence of the experimental 
value of the critical temperature, the Coulomb pseudopo- 
tcntial is unknown. For this reason, we have assumed: 
fi* G (0.1,0.3). The symbol 9 denotes the Heaviside unit 
function and uj c is the cut-off frequency (u> c = 3f2 max ). 



The quantities N(oj) and f(w) denote the Bosc-Einstein 
and Fcrmi-Dirac function, respectively. 

The order parameter function and the wave function 
renormalization factor on the imaginary axis have been 
calculated by using the equations |23l |: 
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The Eliashberg set has been solved for M = 1100. We 
have used the numerical methods presented in the papers 
[26| and [23] . In the considered case, the solutions of the 
Eliashberg equations are stable for T > Tq = 4.64 K. 

In Fig.[TJ the form of the order parameter on the real 
axis has been presented for the selected values of the tem- 
perature and /Lt* = 0.1. Moreover, the rescaled Eliashberg 
function (10a 2 F(fl)) has been also plotted. 

It is easy to sec that for the low frequencies, the non- 
zero values are taken only by the real part of the order 
parameter. The obtained result indicates that in the con- 
sidered range of frequencies the damping effects related 
with Im[A(w)] not exist. Additionally, it has been stated 
that the shapes of the functions Re[A(w)] and Im[A(cj)] 
are correlated with the complicated form of the Eliash- 
berg function. 
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The physical value of the order parameter has been 
calculated by using the equation: 



A (T) = Re [A (w = A(T))]. 



(7) 



In Fig. [2] (A) , the dependence of the order parameter 
on the temperature for the selected values of the Coulomb 
pseudopotcntial has been presented. We notice that the 
value A (T, ) can be parameterized by using the formula: 
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where j3 = 3.4, and: 



A (T ) = 91.93 (p*) — 64.73/x* + 14.47 meV. (9) 

Next, we have determined the dependence of the crit- 
ical temperature on the Coulomb pseudopotcntial. The 
results have been presented in Fig.|2](B). 

It is easy to notice that the value of Tc is high in the 
whole range of the considered values of fj,*. In particular, 
T c G (51.65,20.62) K. 

Additionally, we underline that for large value of the 
Coulomb pseudopotcntial, the critical temperature can 
not be precisely estimated by the classical Allen-Dynes 
or McMillan formula [28|, [2{|. However, the modi- 
fied Allen- Dynes formula derived originally for S1H4 (112)2 
compound works very well [I?} : 
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where /1 and f 2 denote the functions [281 ]: 
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The quantities Ai and A2 take the form: Ai = 2 — 
0.14/**, and A 2 = (0.27+ 10/t*) (V^/W)- 

The second moment of the normalized weight function 
and the logarithmic phonon frequency are given by: 
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The electron-phonon coupling constant has the form: 
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For SiH4 under the pressure at 250 GPa, it has been 
achieved: JZJ2 = 167.64 meV and u>\ n = 72.42 meV. 
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FIG. 2: (A) The dependence of the order parameter on the 
temperature for selected values of the Coulomb pseudopo- 
tential. (B) The critical temperature as a function of the 
Coulomb pseudopotential. The circles represent the results 
obtained by using the Eliashberg equations. The squares and 
triangles are related to the classical Allen-Dynes and McMil- 
lan expression. The solid line has been achieved with help of 
Eq. (10). 
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FIG. 3: The dependence of the ratio Ra on the value of the 
Coulomb pseudopotential. In the inset, the influence of n* on 
A(0) has been presented. 



Now, we consider the low-temperature value of the 
order parameter at the Fermi level (A (0) = A (To)). 
In particular, we have A (0) G (4.56,1.71) meV for 
jU* G (0.1,0.3). The full dependence of A (0) on the 
Coulomb pseudopotcntial has been shown in the inset 
in Fig. El 

Taking into consideration the obtained results, we have 
calculated the dimensionless ratio Ra = 2A (0) /fc^Tc?. 
It has been achieved: Ra G (4.10,3.84) (see also Fig. 
[3]). We notice that the values of Ra for SiH4 compound 
differ significantly from the prediction of the BCS model, 
where R A = 3.53 [H|,[22|. 
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FIG. 4: The dependence of the total normalized density of 
states on the frequency for T — To. The selected values of 
the Coulomb pseudopotential have been assumed. 
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FIG. 5: The dependence of the real and imaginary part of 
the wave function renormalization factor on the frequency for 
T — Tq and T — Tc- In the figure, the rescaled Eliashberg 
function has been also plotted (5a 2 F(Q)). The results have 
been obtained for y* = 0.1. 



The order parameter function on the real axis allows 
to calculate the total normalized density of states f23j : 
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where the symbols N$ and Nn denote the density in the 



superconducting and normal state, respectively. The pair 
breaking parameter T is equal to 0.15 meV. 

The plot of the total normalized density of states has 
been presented in Fig. |4] We have assumed T = Tq and 
the selected values of /i*. It is easy to see the reduction 
of the function N$ (ui) /Nn (ui) by the depairing electron 
correlations. 

In Fig. [5j the real and imaginary part of the wave func- 
tion renormalization factor on the real axis has been pre- 
sented (T = T and T = T c ). 

On the basis of the obtained results, it has been stated 
that Z (w) weakly depends on the temperature in com- 
parison with the order parameter. We can also observe 
the weak correlation between the wave function renor- 
malization factor and the shape of the Eliashberg func- 
tion. 

The dependence of the electron effective mass (m*) 
on the temperature has been determined by using the 
expression: m* = Rc[Z(0)]m e , where m e denotes the 
band electron mass. 

In our case, the results prove that the electron effec- 
tive mass takes the high value in the whole range of the 
superconducting phase's existence. The maximum of m* 
has been observed for T = Tc, where m* = 1.95m e . We 
notice that the Coulomb pseudopotential does not influ- 
ence on the value of [w*] max . 

In the paper, the Eliashberg equations for the SiH4 
compound under the pressure at 250 GPa have been 
solved. We have shown that the value of the criti- 
cal temperature decreases from 51.65 K to 20.62 K if 
H* € (0.1,0.3). The dimensionless ratio Ra takes the 
values beyond the prediction of the BCS model: Ra € 
(4.10,3.84). Additionally, it has been stated that the 
electron effective mass takes the high values and [m*] max 
is equal to 1.95m e for the critical temperature. 
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